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Abstract. A class of exact propagators describing the interaction of an A- level atom 
with a set of on-resonance (5-lasers is obtained by means of the Laplace transform 
method. State-selective mirrors are described in the limit of strong lasers. The ladder, 
V and A configurations for a three-level atom are discussed. For the two level case, the 
transient effects arising as result of the interaction between both a semi-infinite beam 
and a wavepacket with the on-resonance laser are examined. 



PACS numbers: 03.75.Be, 03.75.-b, 31.70.Hq 
The spacetime propagator can be considered as one of the most important tools 
in quantum physics for it governs any dynamical process. However, the knowledge of 
propagators corresponding to non-quadratic Hamiltonians is severely restricted. In this 
line, the spacetime propagator for a (^-potential relevant to tunnelling problems has 
excited much attention [TJ EH HI E] Such interactions turn out to be particularly 
useful to gain physical insight in systems where only integrated quantities are to be 
considered. A thorough discussion of point interactions as solvable models using a 
functional approach can be found in and a formalism to incorporate general point- 



interactions and dealing with different boundary conditions has been developed by 
Grosche (Z1E|- Even though the method is particularly suitable to calculate the energy- 
dependent Green function, a wide class of propagators was derived in such a fashion. The 
incorporation of time-dependent point-interactions has been possible through different 
approaches as Duru's method [0] or the use of integrals of motion (THj. However, most 
of the effort has been focused on the dynamics of structureless particles and to the 
knowledge of the authors no attention has been paid to problems involving internal 
levels. Such state of affairs contrasts dramatically with the current surge of activity in 
atom optics. 

In this paper we use the method of Laplace transform [HE] to tackle particles with 
internal structure. In particular, we shall focus on exact propagators for atom-laser 
interactions, namely, those of an atom interacting with a set of 5-laser on-resonance 
with given interatomic transitions. The method is introduced in section ^ to obtain 
the exact propagator for a two-level atom. Details of the calculations relevant to the 
following sections are here provided. In section Elthe propagators for a ladder, V, and A 
configuration (see Fig. ^) of lasers interacting with a three-level atom are obtained. The 
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general case in which a given state is coupled to an arbitrary number of levels is discussed 
in section |3] where the high intensity limit of the laser is related to state selective 
mirrors. Such kind of systems presents manifold applications in laser coherent control 
techniques such as cold atomic cloud compression JT] , atom mirrors and beam splitters 
[T2] , and different schemes where fast transitions are required as in the implementation 
of logic gates for ion trap quantum computing Moreover, idealised time-of- arrival 
measurements J4j and recently proposed improvements in Ramsey-interferometry with 
ultracold atoms jTHJ rest in a full quantum mechanical treatment of the dynamics of 
such systems. 

1. The two-level atom 

In this section, we use the method of Laplace transform to obtain the propagator 
for a two-level atom incident on a narrow perpendicular on-resonance laser beam. 
Spontaneous decay is assumed to be negligible throughout the paper and we 
shall consider effective one-dimensional systems in which the transverse momentum 
components can be neglected as it is the case for atoms in narrow waveguides |16j . In 
a laser adapted interaction picture, and using the rotating wave approximation, the 
Hamiltonian describing the system is 

where p is the momentum operator conjugate to x, the ground state |1) is in vector- 
component notation and the excited state |2) is (°). The second term in the right 
hand side defines the potential strength matrix V and I2 the two-dimensional identity 
matrix. Equation (Q) may be regarded as the e — > limit of a laser of width e and Rabi 
frequency fl^ keeping Q = Q^e constant. here and in the following is chosen to be 
real. We start then by considering the free propagator for a one-channel problem on a 
Hilbert space of square integrable functions S) (see for instance 

In what follows we shall be interested in describing the dynamics of particles with two 
internal levels. The free propagator (Q = 0) for states on the Hilbert space fj®C 2 is given 
by Ko(x, t\x' , t') = Ko(x,t\x' ,t')l2, in the same interaction picture than Moreover, 
#-type of perturbations can be generally taken into account using the method of Laplace 
transform [Sj which assumes the unperturbed propagator to be known. More precisely, 
the full propagator can be related to the free one through the Lippmann-Schwinger 
equation |H |3] 

K(x,t\x',t') = K (x,t\x',t')-l [ dt" [ dx"Ko{x,t\x",t"W(x",t")K(x",t"\x',t').(3) 

nj t > J -00 

Given that the potential has the form of a point-interaction, the integral over x" 
coordinates is straightforward. One can then take the Laplace transform with respect 
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K(x, s\x', 0) = K (x, s\x', 0) - ^ ( ~ ° KqX,s£, 0) 

ft \ A (x,s|f,0) 



K(e,s|x',o), 



(4) 



where we have made use of the convolution theorem (£[/ * g(t — t')f{t')dt'\ = g(s)f(s)). 
By setting x = £ it is explicitly found that 

-l 



1 



f^o(0,s|00) 



f^ (0, S |00) 



K o (e,s|x',0). 



(5) 



Next, we note the exact expression for the Laplace transform of the single-channel free 
propagator (j2J, 



K (x, s\x', 0) 



m 



lifts 



(6) 



which becomes necessary for evaluating the inverse of the matrix. Combining (j3J) and 
(J5J) the Laplace transform of the exact full propagator is obtained, and taking the inverse 
transform one can find the spacetime propagator 



K(x,t\x',0) = K (x,t\x',0) - 
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Fortunately in our case, the resulting matrix element can be related after taking partial 
fractions with standard results so that 

mV 



K(x, t\x, t ) = K (x, t\x', t') - 



Ah 2 



EamVo(\x-£.\ + \£-x' \) ,-mVgt 
e h 2 2ft 3 " 

a=± 



x erfc 



a 



imVnt 1 2m .. .. .. ... 



a 1 
1 a 



(9) 



A more compact expression can be obtained rewritting the full propagator in terms of 
the Moshinsky function (see | Appendix A] ), 



K(x, t\x ,t) = K (o:, tlx , £ 
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with k = —imVo/h 2 for short. One should notice that this propagator opens up the way 
to a whole variety of problems involving quantum dynamics of two-level atoms. A most 
relevant fact is that as long as the full Hamiltonian can be written as a direct sum, the 
same property holds for the propagator. Therefore, if H = © s=1 H s , 



K(x, t\x', t') = Q)K s (x,t\x',t'), 



:iv, 
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Figure 1. Different configurations of a three-level atom interacting with a pair of 
on-resonance fields: a) ladder, b) V, and c) A type. 



which becomes very useful when one wishes to study the dynamics in a given subspace. 
In sections 0] and we shall focus on some analytical examples dealing with the quantum 
dynamics of semi-infinite beams and wavepackets of 2-level atoms. 

2. The three level atom 



The procedure to obtain the propagator can be extended to atoms of many-levels. Even 
though the general expression can be found, the calculation is straightforward up to the 
inverse Laplace transforms. In this section we consider a three-level atom subjected to 
two on-resonance lasers in different configurations (see Fig. EJ. First, we shall focus 
on the ladder configuration which involves the transitions |1) — > |2) and |2) — > |3) by 
means of lasers located at the same position. The Hamiltonian describing the system is 
then 

/ Q 12 \ 



H c = JLl a + VS(x - = f-l 8 + -S(x - 
2m 2m 2 



Q 12 Q 23 

n 23 



.(12) 
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Following section ^one finds 
K(x,t\x',t') = K (x,t\x',t') - -^/T 1 

h 



K V(x, 0) ( 1 3 + -K o V(0, s\0, 0)) K (£, s\x', 0) 



Plugging the explicit form of the free propagator evaluated at the different positions, 
taking partial fractions, working out the inverse Laplace transform and rewritting the 
result in terms of the Moshinsky function, one finds the exact expression for the full 
propagator to be 



K(x, t\x', t') = K (x, t\x>, t f )l 3 - ^ J2 M ( \ x ~ *l + K - ^ -^ ^/^ ^ 
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Another relevant case is the V configuration in which a pair of lasers couples the 
state |1) with both \2) and |3) levels, in such a way that the resulting Hamiltonian is 
given by 



/ 
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(14) 



The exact propagator can be found to be 
K(x 7 t\x\t') = K (x,t\x',t')l 3 -^J2 M 
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From the mathematical point of view such configuration is actually very similar to A 
scheme described by the Hamiltonian 
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whose propagator reads 
K(x,t\x',t') = K (x,t\x',t')l 3 
( 
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One may then study the space-time dynamics of coherent trapping which results from 
the destructive quantum interferences between the two transitions, for initial states in 
a superposition of the two lower levels |1) and |2) [T5] . 



3. The iV-level atom 

In this section we extend the previous approach to an N- level system living on <S> C N . 
Suppose that an iV-level atom is subjected to the action of N on-resonance lasers 
{Mlij/2 = Vi) all of which are located at the same position and couple the \j) level 
(1 < j < N) with the (N - 1) levels as shown in Fig. ©, 
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Figure 2. Configuration of a iV-level atom interacting with a (N — 1) on-resonance 
fields coupling the |j) to any other level. 
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Simplified configurations can be obtained just by setting certain coupling elements to 
zero. The exact propagator for such system reads 

K.(x, t\x', t') — Ko(x, t\x', t')ljv 



2h 2 



M ( \x — £| + |£ — x'\, —iam^, ^ 
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V N V m . (19) 
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aV 2 ) 

where V m = \J^^vf ■ Incidentally, the ladder configuration in which only successive 
levels are coupled to each other, does not admit a similar generalisation to iV-level given 
the complexity of the inverse matrix. 



aV N Vj-i V N V m aV N V j+ i 



3.1. Intense lasers and propagators for state-selective mirrors 

For the one-channel case, the effect of an infinite strength of the ^-potential is 
tantamount to a hard- wall boundary condition This feature has been extensively 

exploited in the exact perturbation theory developed by Grosche to obtain a wide class 
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of Green functions and propagators including boundaries (say, in half-space or boxes) 
[3IH]. In atom optics, it is also well-known that a very intense laser behaves as a totally 
reflective mirror for the states coupled by it. Therefore, we shall consider the limit in 
which one of the lasers is made infinitely strong. In order to do so, it is convenient to 
note that the following relation holds [5 , 




hm — £ — £ + £ — x\, —iam— — — ,- =aK (x - £ + £ - x ' ,t 



0,0 



The general expression for the propagators with iV-lasers in which the n-th one (coupling 
the states |j) and |n}) is made infinitely strong becomes diagonal, 

[K(x, t\x', t% k = K (x, t\x', t')8 ik - K (\x - £|, 1 1 W - £|, 5 ik (5 jk + 8 nk ), (20) 

5ij being the Kronecker delta. The laser behaves then as a totally reflecting mirror 
selective to the sates coupled by it, and suppresses for such states any possible excitation 
due to the presence of other lasers located at the same position. 

Another interesting case is the on-resonance excitation of a couple of levels with 
high-intensity lasers. This leads to the limit where the strength coefficients of the two 
lasers go to infinite, Vi,V n — > oo and n > I, but its ratio is kept constant, Vi/V n = c 
with c > 0. For such case one finds 



K(x,t\x',t% k = -K (\x - £\,t \x' - £\,t'j -jJ= = (5 ni 8i k + 5 u 8 nk ) 



+ 5: 



ik 



K (x,t\x',t') - K (\x -£\,t \x' - £|,/) (^5 kj + -j^==5 ik + -j2==5 nk ^j 



.(21) 



The difference now arises from the fact that the state-selective mirrors have a finite 
reflectivity for states |Z) and \n), to which excitation is allowed. However, notice that 
for the \j) state the laser still mimics a totally reflecting mirror. 



4. Moshinsky shutter 

We next study a time-dependent multi-channel scattering problem and consider a 
monochromatic beam of two-level atoms in its ground-state incident on a totally 
absorbing shutter which is suddenly removed at time t = 0. Such kind of setup is 
usually referred to as a Moshinsky shutter, ever since the seminal paper [201 which led 
to the discovery of diffraction in time. The conditions on the reflectivity of the shutter 
can be easily modified to more general cases EI] The initial state is then of the 
form 

V(x,t = 0) = e ikx Q{-x)\l), (22) 

where G(x) is the Heaviside step function. Equation ()22j) is an obvious generalisation 
of the Moshinsky type of initial condition for a single-channel problem, discussed in the 
context of diffraction in time. We note that such kind of state is not normalisable in 
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Figure 3. Probability density of the ground and excited states, total probability 
density in the presence of the laser and for the free case, 50 ms after removing the 
shutter initially located at 50 fim. The incident beam moves at v = 1 cm/s (mass 
of 87 Rb). The picture is taken at the instant in which the classical profile (the step 
function <d(vt — x)) reach the laser. 



the usual sense, yet accurately describes certain experimental setups [22] and provides 
a basis for wavepacket analysis. 

The time evolution can be studied using the superposition principle. If we consider 
first the free evolution, with Q = 0, then the result is that of diffraction in time, 

/oo 
dx'K (x,t\x',t')^f{x,t = 0) 
-oo 

= M(x,k,ht/m)\l). (23) 

This solution was found by Moshinsky and has been observed in a wide variety of 
experiments with ultracold neutron interferometry for the one-channel case [221 ■ It is 
a well-known fact that it tends with increasing time to the stationary wavefunction. 
We further notice that due to the absence of coupling between the internal states in 
K. (x,t\x',t'), the excited state |2) is not populated, in agreement with (fTTf . 

An alternative configuration in which the beam tunnels through a 5-barrier was 
discussed in jU EHl EH • 

Let us now look at the time evolution in the presence of the laser. Using the integral 
u4]) in Appendix A[ the exact solution can be found in close-form, 



/oo 
dx'K(x,t\x',t')&(x,t = 0) 
-oo 

=*•<*■«)+ 5 1: Q 

m V- * ' 



a=± 

x \ [M(\x - £| + £, k, ht/m) - M(\x - £| + £, a«, M/m)]}.(24) 
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Figure 4. Total probability density for an incident beam in the ground state with 
the laser and for free space, 150 ms after removing the shutter initially located at 100 
(im. The incident beam moves at v = 1 cm/s (mass of s7 Rb), exhibiting interference 
for all x < £, as a result of the reflexion from the laser. 
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Figure 5. Space-time density plot of the population in the excited state, |(2|'4')| 2 , 
with v = 1 cm/s, and the position of the laser £ = 200 /im. The grey scale changes 
from dark to light as the function values increase. 

The total probability density is plotted in Fig. El The velocity of the incident beam 
is chosen in all simulations to satisfy v = hq/m = Vq/H, for which one finds, after solving 
the two-channel stationary problem, that the reflexion and transmission probabilities in 
the excited state, I-R2I 2 = 1^2 1 2 , are maximised and indeed equal those in the ground 
state, |-Ri| 2 = I T\ | 2 = 1/4. The paradigmatic oscillations on the probability density, 
main feature of the diffraction in time described by the Moshinsky function (|^o| 2 ), 
is modified for all x < £ due to the interference which arises with the reflected part. 
Indeed, for later times such interference completely dominates and the density profile is 
dramatically perturbed as shown in Fig. HJ 
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In Fig. [5] we plot the time evolution of the probability density in the excited state 
1 2). The delta laser is shown to behave as a point-like source of atoms. Moreover, the 
pattern exhibits diffraction in both time and space domain. 



5. Wavepackets dynamics 

Many experiments deal with finite samples rather than beams. In the one-channel case, 
the tunnelling dynamics of wavepackets through narrow barriers have been examined 
in a series of works pEJEJlEE]. In addition, the phenomenon of quantum deflection was 
predicted in the presence of semi-transparent and perfect mirrors [2l)] |2*7| l2*%] . 

We next consider normalisable states belonging to £ 2 (M) <S> C 2 . In particular we 
study the dynamics of eigenstates of a hard-wall trap which are released at time t — 
and launched with momentum hq against the on-resonance delta laser. All-optical box 
traps have been recently been obtained in the laboratory j2H] • More precisely we assume 
an initial sine-wavepacket given by 

*om = o) = ^a/Ie^^i 1 ) ( 25 ) 

with q n p = q + (3nn/L. The expansion of such state has recently been discussed at 
the single-channel level in free space j3U], i n the presence of gravity [51], and when 
generalised to the many-body Tonks-Girardeau regime [32]. For the interaction with 
the on-resonance delta laser one can actually propagate in time this initial condition. 
The time evolved wavefunction is 

/oo 
dx'K Q {x,t\x',t = 0)*(x',i' = 0) 
-oo 

= ^^^P[e^ L M(x-L,q nP ,ht/m) - M(x,q n(3 ,M/m)]\l) (26) 



p=± 



and, 



a,(3=± V 7 



k — an 

x h ip e L/h [M(\x - f | + £ - L, k, ht/m) - M(\x - £| + £ — L, an, ht/m)} 

_ [ M (\x - ei + e, fe, ^M) - m(\x - ei + e, ^/^)] }■ ( 2 7) 



Figure [H] shows the dynamics when the incident momentum is such that maximum 
excitation is achieved. The laser acts as a beam splitter dividing the wavepacket into 
two parts, the transmitted one being similar in shape to the freely evolving wavepacket, 
whereas the reflected part exhibits interference. Notice that ([2%]) and ([2*6^) admit a 
simple generalisation for the respective iV-level problem. In order to do so it suffices to 
consider the suitable momenta and prefactor of the matrix in the propagator for each 
of the channels. 
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Figure 6. Ground, excited and total probability densities for an incident sine- 
wavepacket (n = 1) released at t = from a hard wall trap of size L = 50/im and 
centered at 25/im from the origin, with q — mVo/H 2 and velocity lcm/s, after 100 ms 
of evolution. The position of the laser is £ = 100/xm. The freely time evolved sine 
wavepacket is also plotted. 



6. Discussion 

We have generalised the method of Laplace transform to include point-like perturbations 
in the dynamics of particles with internal structure. In such a fashion we have obtained 
the spacetime propagator for an A-level atom interacting with a set of on-resonance 
delta lasers. For strong lasers, the propagators for state-selective mirrors have been 
obtained. A similar procedure could be applied to the Green function for which an 
exact perturbation theory has been developed and extensively discussed by Grosche |H] . 
In the 3-level atom case, the ladder, V, and A configuration have been considered. The 
inclusion of time dependence in the propagators to adiabatically turn the lasers on and 
off independently from each other is an open problem which would provide access to 
the space-time dynamics of coherent population trapping and other relevant phenomena 
[TH] . The dynamics of a semi-infinite beam and a wavepacket of two-level atoms incident 
on the laser, with a straightforward generalisation for the A-level case, have also been 
worked out. 
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Appendix A. The Moshinsky function 

The Moshinsky function arises in most of the problems where ID quantum dynamics 
involves sharp boundaries well in time or space domains. Similarly it is found when 
considering free propagators perturbed with point-interactions. 
Its standard definition reads 

if 

M(x, k, t) := ^w(-z), (A.l) 



' ' r (*--), (A.2) 



where 

1 

z = - 

2 v V t 

and the so called Faddeyeva function [3H] w is explicitly defined as 

_ 2 If e"" 2 

w(z) := e z erfc(-iz) = — / du , (A. 3) 

ITT J r u — z 

T_ being a contour in the complex z-plane which goes from — oo to oo passing below 
the pole. After [20J, M(x, k, r) has been named the Moshinsky function. 

For the exact time evolution in the presence of a laser with a cut-off plane wave or 
hard-wall eigenstates as initial conditions, we find integrals of the form 

/x ikx 
alx'e ikx 'M{ax' + b, c, r) = — \M{ax + b, c, r) - M{ax + b, -k/a, r)l . (A.4) 

i(k + ca) 

For more details we refer the reader to ^j. 
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